Pseudo-Rabi oscillations in superconducting flux qubits in the classical regime 
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Nonlinear effects in mesoscopic devices can have both quantum and classical origins. We show 
that a three- Josephson-junction (3JJ) flux qubit in the classical regime can produce low-frequency 
oscillations in the presence of an external field in resonance with the (high-frequency) harmonic 
mode of the system, lu. Like in the case of quantum Rabi oscillations, the frequency of these pseudo- 
Rabi oscillations is much smaller than u> and scales approximately linearly with the amplitude of 
the external field. This classical effect can be reliably distinguished from its quantum counterpart 
because it can be produced by the external perturbation not only at the resonance frequency ui and 
its subharmonics (u)/n), but also at its overtones, nui. 



The many advances in understanding the quantum behaviour of meso- and macroscopic systems over the last decade 
are due, to a large extent, to the investigation of superconducting qubits in the context of quantum computing [H, @]. 
These led to the demonstration of their quantum behavior, at the level of one to four qubits [!, Bl, H, IE 0] • 

Since flux qubits can exhibit quantum superpositions of states differing by a macroscopic number of single-electron 
states, and the relevant observables can be easily accessed experimentally, these devices provide unique opportunities 
to investigate the quantum-classical frontier. Further scaling up of superconducting qubit networks is thus motivated 
by the needs of quantum information processing in solid state as well as by the fundamental interest in probing the 
limits of the applicability of quantum mechanics. 

Once measurements are taken into account, the quantum behaviour becomes essentially nonlinear. However, already 
at the classical level, nonlinearities are unavoidable in superconducting qubits, due to the nonlinear behaviour of 
Josephson junctions. It was recently pointed out that in a phase qubit, which is a biased single Josephson junction, 
the classical nonlinearity can produce effects with characteristics similar to Rabi oscillations [8j and Ramsey fringes 
, which are often considered as signatures of quantum behaviour in two- level systems 

The results of [H, Q do not undermine the common understanding that superconducting qubits demonstrate 
quantum behaviour, since the latter was established by a set of independent methods Rather it attracts our 
attention to the interesting possibility of coexistence of similar nonlinear classical and quantum effects in the same 
device. How to distinguish classical versus quantum behavior in these qubits is an important question and the focus 
of this work. 



Rabi oscillations (e.g., Ref. llj, p. 89) are coherent quantum transitions in a two- level system driven by an external 



ac field of amplitude A and with the characteristic frequency 

n = VA 2 + (M 2 , (i) 

where 6u> = \lu — luq\ is its detuning from the interlevel distance, ojq. In resonance, f2 — A. This linear dependence 
on the field amplitude, and Q being much less than other characteristic frequencies in the system, help identify Rabi 
oscillations. Indeed, they were observed in all types of superconducting qubits 0, and their appearance is sometimes 
used as decisive evidence in favour of quantum behaviour. Multiphoton Rabi oscillations, at loq — nuj, were also 
observed [12J • 

In this paper we investigate a three- Josephson-junction (3JJ) flux qubit in the classical regime. With two inde- 
pendent variables instead of one, this is a richer system, than the phase qubit of [8j. We find that the resonant 
high-frequency driving uj produces low-frequency oscillations of the magnetic flux which are very similar to Rabi oscil- 
lations. We also show that a qualitative difference exists between these two effects, which allows to reliably distinguish 
them in experiments. 

Following Ref. [l3j |. we consider a 3JJ flux qubit in the limit of negligible self-inductance L — * 0. The magnetic flux 
through the loop then equals the external applied flux <I> e , and we introduce the reduced flux ip e — 27r$ e /$ - Here 
$o = h/2e is the flux quantum. Bistability is achieved due to the presence of three junctions in the loop. Due to 
the single-valuedness of the superconducting wave function, one of the phase shifts across the Josephson junctions is 
eliminated through <pi + (f2 + ^3 = ¥>ei leaving two independent variables, = (epi + </?2)/2 and \ = ifi ~~ ( /'2)/2. 



2 




In the classical regime, the phase dynamics of the ith Josephson junction (i = 1, 2, 3) can be described by the RSJJ 
model 14] , in which the current is given by 



T hC l d 2 



h d 

2eRi ~dt 



ipi + I cl sin ipi 



(2) 



Here Cj is the junction capacitance, Ri its normal resistance, and I C i its critical current. We now neglect the effect 
of thermal fluctuations, which are not crucial here. Following the common choice of parameters, we set C\ — C% — 
C, I c i = I C 2 = la Ri = R"2 = R', C3 = otC, I C 3 = al c , i? 3 = R/a; (0.5 < a < 1). Introducing the dimensionless units, 
LUot — t, 7 — hujo/2eRI c = ll>o/lur, where loq = y (2e/ c )/(ftC), and ujr — (2eRI c )/h, we can write the equations of 
motion for the variables 9, x : 



-r^X + 7-i-X = - cos sin x, 

(l + 2a)-^0 + 7 (l + 2a)-^-0 = 

d 2 d 
sin 9 cos x + a siri(ip e — 29) + a-—r(p e + aj— — ip e . 

dr z ar 



(3) 



We consider the external flux 

The energy of the system is thus 
H = Ej 



ip e (t) = 27r($ e /$ ) = ft + ^sin(wr). 



cos cos x — ^cos((p e — 29) 



(4) 



(5) 



where Ej = HI c /2e is the Josephson energy. The canonical momenta are p x = Ejd\/dr, 
effective potential is given by (Fig. [1]) 



Ej(l + 2a)d9/dT. The 



U(9, x) = - cos cos x ~ a cos(<p e - 29) /2. (6) 

If the dc static bias is ipf = 7T, the system has degenerate minima at O = ± arccos [(2a) -1 ] , xo — 0- The "plasma" 
frequencies of small oscillations around them are (in units of ujq) ujb = yfl — (2a) _1 , uj x — (2a)~ 1 / 2 . 

In the presence of the external field (0} the system will undergo forced oscillations around one of the equlibria. For 
a = 0.8 [IH, which is close to the parameters of the actual devices 0, 0, EH , the values of the dimensionless frequencies 
become lob ~ 0.612, and lo x w 0.791. Solving the equations of motion ([3]) numerically, we see the appearance of slow 
oscillations of the amplitude and energy superimposed on the fast forced oscillations (Fig. [5]), similar to the classical 
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FIG. 2: (a) Driven oscillations around a minimum of the potential profile of Fig. [T]as a function of time. The driving amplitude 
is ip% — 0.01, driving frequency to = 0.612, and the decay rate 7 = 10 -3 . Low-frequency classical beat oscillations are clearly 
seen, (b) Low-frequency oscillations of the persistent current in the 3JJ loop, (c) Same for the energy of the system. 



oscillations in a phase qubit (Fig. 2 in [8|). The dependence of the frequency of these oscillations on the driving 
amplitude shows an almost linear behaviour (Fig. [3]), which justifies the "Pseudo-Rabi" moniker. 

A quantitative difference between this effect and true Rabi oscillations is in the different scale of the resonance 
frequency. To induce Rabi oscillations between the lowest quantum levels in the potential ^ , one must apply a signal 
in resonance with their tunneling splitting, which is exponentially smaller than ujq. Still, this is not a very reliable 
signature of the effect, since the classical effect can also be excited by subharmonics, ~ o-io/ n , as we can see in Fig. |4] 

The key observable difference between the classical and quantum cases, which would allow to reliably distinguish 
between them, is that the classical effect can also be produced by driving the system at the overtones of the resonance 
signal, ~ nujQ (Fig. 0]). This effect can be detected using a standard technique for RF SQUIDs 17]. The current 
circulating in the qubit circuit produces a magnetic moment, which is measured by the inductively coupled high-quality 
tank circuit. For the tank voltage Vr we have 



d 2 V T 
dt 2 



1 dVr 9 T r 2 » r dlq 
tt dt dt 



(7) 



where tt = RtCt is the RC-constant of the tank, lot = (LtCt) its resonant frequency, M the mutual inductance 
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FIG. 3: The dependence of the pseudo-Rabi frequency on the driving amplitude if>% for to — 0.6, 7 = 10 3 . The solid line, 



^ = 0.35 [(^) 2 + (w-0.63) 



21 1/2 



is the best fit to the calculated data. 




FIG. 4: The average energy H of the system as a function of the driving frequency u). The main peak (u>o sa 0.6) corresponds 
to the resonance. The left peak at cjo/2 is the nonlinear effect of the excitation by a subharmonic, similar to a multiphoton 
process in the quantum case. The right peak at 2ojo is the first overtone and it has no quantum counterpart. Here (p^ = ir; ip" = 
0.05; 7 = 10" 3 . 
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between the tank and the qubit, and I q (t) the current circulating in the qubit. The persistent current in the 3JJ 
loop can be determined directly from Its behaviour in the presence of an external RF field is shown in Fig. [2};. 
Note that the sign of the current does not change, which is due to the fact that the oscillations take place inside 
one potential well (solid arrow in Fig. [T]), and not between two separate nearby potential minima like in the quantum 
case. (Alternatively, this would also allow to distinguish between the classical and quantum effects by measuring the 
magnetization with a DC SQUID.) 

There can be two strategies in detecting the slow oscillations, as long as the measurement time is smaller than their 
decay time. First, one can directly measure the time-dependent voltage in the tank circuit, which from ([7]) is 

Vt{uj) = — ^ I q (u). 8 

CO 2 - Uy + ILU I T T 

Alternatively, one can measure the spectral density of the signal in the tank, 

a uPu^M* 2 

{Vt) "~ (^-c4 )2 + cj2/t 2 1V- i J ) 

Choosing the tank frequency lot close to the classical low "Rabi" frequency f2, in either case we use the tank as 
a filter, which removes the interference from the large, high-frequency driving field. Note that in the first case we 
have an additional discriminant of the quantum versus classical behaviour of the system. In the classical case, the 
time-dependent regular oscillations of the tank voltage or current due to the oscillations in the qubit, ((Sj), can be 
measured directly. In the quantum case, this is impossible, due to the uncertainty principle, and the oscillations are 
only seen in the statistics of the measurements. Remarkably, the effect on the correlators can be directly observed in 
both quantum and classical cases (in the quantum case such measurements are not limited by the decay time jl6|). 

In conclusion, we predict that a 3JJ flux qubit driven by a resonant external field will exhibit classical low-frequency 
oscillations superficially-similar to the quantum Rabi oscillations in a driven two-level system. Both effects can coexist 
in the same region of parameters. A qualitative difference between the two (allowing to reliably distinguish between 
them) is that the classical effect can be driven also at the overtones of the resonant frequency. 
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